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ABSTRACT 



We report IRAM millimeter interferometry of three z ~2.4 to 3.4, SCUBA 
deep field galaxies. Our CO line observations confirm the rest-frame UV/optical 
redshifts, thus more than doubling the number of confirmed, published redshifts 
of the faint submm-population and proving their high-z nature. In all three 
sources our measurements of the intrinsic gas and dynamical mass are large (10^° 
to 10^^ M0). In at least two cases the data show that the submm sources are 
part of an interacting system. Together with recent information gathered in the 
X-ray, optical and radio bands our observations support the interpretation that 
the submm-population, at least the radio detected ones, consists of gas rich (gas 
to dynamical mass ratio ~0.5) and massive, interacting starburst/AGN systems. 

Subject headings: cosmology: observations - galaxies: formation -galaxies: high- 
redshift - galaxies: evolution 
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1. Introduction 

The extragalactic far-IR/submm background is probably dominated by luminous in- 
frared galaxies (LIRGs/ULIRGs: Ljr ~ 10^^-^ to 10^^ L©) at z > 1 (e.g. Small et al. 1997; 
Bertoldi et al. 2002; Scott ct al. 2002; Cowie et al. 2002). However, far-IR/submm sources 
in most cases have relatively poorly known positions and frequently have only weak counter- 
parts in the rest- frame UV and optical (Small et al. 2000, 2002; Dannerbauer et al. 2002). 
As a result redshifts have thus far been confirmed with CO interferometry for only two of the 
~ 100 detected systems (Blain et al. 2002). Recently a subgroup of the authors have obtained 
optical spectroscopic redshifts for a number of sources detected with the SCUBA camera at 
850;um, to a large extent aided by more precise positions derived from deep 1.4 GHz VLA 
observations of the same fields. Here we report the first results on the millimeter CO line 
follow-up of these submm sources. We believe these observations mark a sensitive break- 
through in the notoriously difficult study of the faint far-IR/submm galaxy population. 



2. Observations 

The observations were carried out between late summer 2002 and winter 2003 with 
the IRAM Plateau de Burc interferometer, consisting of six 15m-diameter telescopes. We 
used the compact D configuration, and for follow-up observations of two of the sources the 
more extended EC configurations. The correlator was configured for CO fine and continuum 
observations to simultaneously cover 580 MHz in the 3 mm and 1.3 mm bands. The frequency 
settings were adjusted for all three sources to optimize the CO line centering in the bandpass. 
SMMJ04431+0210 was observed between Sep '02 and Feb '03 in D and BC configurations 
for a total integration time of 22 hrs. SMMJ0943 1+4700 was observed in D configuration 
only, in Nov '02, for 13 hrs. SMMJ16368+4057 was observed between Sep '02 and Feb '03 for 
24 hrs. All sources were observed in very good observing conditions. We calibrated the data 
using the CLIC program in the IRAM GILDAS package. Passband calibration used one or 
more bright quasars. Phase and amplitude variations within each track were calibrated out 
by interleaving reference observations of nearby quasars every 20 minutes. The overall fiux 
scale for each epoch was set on MWC 349. After fiagging bad and high phase noise data 
(less than 2% at 3 mm), we created data cubes in natural weighting. 
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3. Results 

Figures 1 and 2 show the CO spectra and maps. The derived properties are hsted in Ta- 
bles 1 and 2. We adopt a flat, Qu = 0.3 A-cosmology with Hq=70 kms~^Mpc~^. To convert 
CO luminosities to gas masses, including a 37% correction for helium, we adopt, under the 
assumption of constant brightness temperature for the lowest rotational transitions from 1-0 
to 4-3, a factor a—O.S M0/(K kms~^ pc^)=0.2q; (Galactic), as derived from observations of 
z ~ 0.1 ULIRGs (Downes & Solomon 1998). The gas masses are probably uncertain by a fac- 
tor of at least 2. We estimate dynamical masses from Mdyn sin^i (Mq) = 4 x IO^Aw^^hm^- 
Here we assume that the gas emission comes from a rotating disk of outer radius R (kpc) 
observed at inclination angle In a merger model the dynamical masses would be a factor 
of 2 larger (Genzel et al. 2003). The numerical constant incorporates a factor of 2.4 between 
observed FWHM velocity width of the line emission, AfpwHM (km/s), and the product of 
rotation velocity and sim. This factor is estimated from model disks taking into account local 
line broadening, beam and spectral smearing. We deduce IR luminosities from the SBO/im 
continuum flux densities S by adopting a modifled greybody model (T=40 K) with oc i/^-^ 
dependent emissivity, such that in the range from z—2 to 3.5 Lir, (L©) = 1.9 x 10^^ Sgso with 
Sgso in mJy (Blain et al. 2002). These luminosities are uncertain by a factor of 2 to 3 since 
dust temperature and emissivity law may vary from source to source (Blain et al. 2003). 
In the following, all linear sizes, masses and luminosities are corrected for the foreground 
lensing factors in Table 2. 

SMMJO443I+O2IO (^=2.51) was originally found in the SCUBA Lens Survey (S85o=7.2 mJy; 
Smafl et al. 1997; 2002). It is located behind the ^=0.18 cluster MS0440+02. Small et al. 
(1999) identifled the submm source with the K = 19 A extremely red object (ERO) N4 about 
3" NW of an edge-on cluster spiral galaxy Nl (N4: R-K=6.3; Prayer et al. 2003). Prayer 
et al. deduced a redshift of z=2.5092±0.0008 from Ha/[Nll]/[0 ill] line emission. The rest 
frame optical line ratios suggest that N4 is a composite star burst/narrow line AGN. We 
adopt the foreground lens magniflcation of 4.4 deduced by Small et al. (1999). Our BCD 
conflguration data show a strong CO 3-2 line centered at 2;=2.5094±0.0002, -|-17 (± 17) 
kms~^ redward of the nominal redshift of the Ha line. The line width is PWHM 350 ± 
60 kms~^, somewhat smaller than that of Ha (520 kms~^). The integrated CO line flux 
corresponds to a total gas mass of 8x 10^ M0. Most of the CO line emission comes from within 
1". The CO emission centroid is 1.1" SW of the near-IR position of N4, as determined by a 
new astrometric solution of the near-IR/radio astrometry we obtained by comparing USNO 
stars with radio sources in the field (N4: 04*^43^07.25' 02°10'24.4" (J2000); the uncertainty 
is ±0.5"). This new position of N4 is 2.2" E and 0.6" S of the position reported by Small 
et al. (1999) based on the APM coordinate system. The 1.3 mm continuum data show a 
marginally significant detection (1.1 mJy, 3.7a) near the position of N4. We set a limit to 
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the CO 7-6 emission of < 0.8 Jy kms~^ (2(j)- For comparison, the Ha emission exhibits 
a velocity gradient of > 400 kms~^ over about 1" and along the slit at p. a. —14°, that is, 
at about 50° relative to the direction of the extended CO emission. The centroid of the 
Ha emission is on N4. It thus appears that the rest-frame submm and optical observations 
sample a similar region (size ~1") but with some differences in the spatial structure in the 
two wavelength ranges. SMMJ04431+0210 has by far the lowest intrinsic IR luminosity 
(3 X 10^^ L0) and gas/dynamic mass of our three galaxies. We deduce an upper limit to the 
dynamical mass of 4.5x10^ sin~^ i Mq for a source diameter of < 1". In terms of luminosity, 
gas and dynamical mass, SMMJ04431-I-0210 thus resembles local ULIRGs. 

SMMJO943I+47OO (^=3.35) was first identified by Cowie et al. (2002) in a deep SCUBA 
map of the 2;=0.41 cluster Abell 851 (S85o=10.5 mJy). They estimated a foreground lens 
magnification of 1.2. Ledlow et al. (2002) proposed that the counterpart of the SCUBA 
source is the 1.4 GHz radio source H6 (72/iJy), for which they identified a redshift of 2;=3.349 
from Lya. H6 appears to be a UV bright, narrow line Seyfert 1 galaxy. We find strong CO 
4-3 emission with a flat-topped profile and FWHM 420 kms~^ centered at z=3.3460 ± 
0.0001, —207 (±7) kms^^ blueward of the nominal Lya redshift. The 1.3mm continuum 
was also detected with 2.3 ± 0.4 mJy and is centered at the same position. CO line and 
continuum emission are centered 3.8" W and 1" S of the position of H6 (24 kpc in the 
source plane), positionally coincident with the second, weaker 1.4 GHz source H7 (55 /xJy) 
at 09*^43^03.7' 47°00'15.1" (J2000) (Fig.2). H6 and H7 are very probably physically related. 
The gas mass deduced from the CO 4-3 fiux is 2.1x10^° Mq, and the IR luminosity based 
on S850 is 1.7x10^^ Lq. For an assumed source size of 1" of H7 we infer a dynamical mass of 
2.5xl0^°sin"^(i) M©. A lower limit to the virial mass of the H6/7 system is 6x10^° Mq. 

SMMJ16368+4057 (z=2.39, ElaisN2 850.4) was identified by Ivison et al. (2002; Sgso^ 
8.2 mJy) in the 8mJy SCUBA blank field survey of the Elais N2 field (Scott et al. 2002) 
with a 220/xJy bright 1.4 GHz radio source. Optical spectroscopy by Chapman et al. (2003) 
and Small et al. (2003) showed bright Lya, Nv, Civ, [Oil] and [Olll] emission with a 
complex spatial and velocity structure. Small et al. (2003) proposed that N2 850.4 consists 
of a UV bright starburst galaxy at 2;=2.380 (±0.002), plus a Seyfert 2 galaxy at 2;=2.384 
(±0.003). There is no evidence for gravitational lensing. We also detected the 1.3mm con- 
tinuum emission (2.5 ± 0.4 mJy) and CO 7-6 emission. The CO 3-2 emission is very broad 
(840 kms-^ FWHM) and is centered at ^=2.3853 (±0.0004), ±115 kms"^ (±36 kms-^) red- 
ward of the nominal redshift of the Sey2 nucleus. CO 7-6 emission is tentatively detected in 
a narrow component centered — 200kms~^ of the CO 3-2 line centroid. At the 7-6 peak the 
observed 7-6/3-2 brightness temperature ratio is 0.5 ± 0.2. Large velocity gradient modeling 
of this ratio indicates that the higher excitation CO 7-6 emission may come from a specific 
warm (T > 50 K) and dense (n(H2) > 10^ cm~^) region. The absolute astrometry of the rest- 
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frame optical/UV, submm and radio positions (each ±0.3") is not yet sufficient to establisfi 
with certainty the relative locations of the emission sources at different wavebands. Relative 
positions are more precise and indicate that the UV bright source is about 0.5" W or SW of 
the optical source, and both have a size of about 0.7", or 5.7kpc. Likewise, we find that the 
different submm components are spread over ~ 0.7", with the CO 7-6 and submm continuum 
about 0.3" to the SW, while the CO 3-2 emission is centered to the NE. Keeping in mind 
that these differences are marginally significant we note that the spatial offsets appear to 
be along p. a. 45°, the direction of the separation between the UV and optical hne emission 
sources (Small et al. 2003). The 1.3 mm continuum (1.6 x 10^^ L©) is unresolved, with an 
upper hmit of about 1". The gas mass estimated from the CO emission is about 5.4 x 10^° 
Mq, and the dynamical mass is 4.5xl0^°sin~^i M© for an adopted source diameter of 0.7". 

4. Discussion 

Our observations of three SCUBA selected galaxies confirm the redshifts identified from 
rest-frame UV/optical spectroscopy, although for SMMJ09431-I-4700 the optical redshift is 
inferred from a source that is physically distinct from the submm source. Our data more 
than double the number of published, mm-confirmed SCUBA redshifts. In addition to the 
three galaxies discussed here, these are SMMJ14011+0252 (2;=2.56; Fraycr ct al. 1999; Ivison 
et al. 2001; Downes & Solomon 2003) and SMMJ02399-0136 (^=2.81; Frayer et al. 1998; 
Ivison et al. 1998; Genzel et al. 2003). Our observations confirm that at least some SCUBA 
galaxies are luminous and gas-rich systems seen at a similar epoch to the UV-bright QSO 
and Ly-break galaxies populations (Boyle et al. 2000; Steidel et al. 1999). 

All five SCUBA galaxies are rich in molecular gas. For the CO luminosity to gas mass 
conversion factor appropriate for local ULIRGs the median gas mass of the five SCUBA 
sources is 2.1 (±1.7) x 10^° M0, similar to the median molecular gas masses found in high-^ 
QSOs (e.g. Alloin et al. 1997; Downes et al. 1999; Guilloteau et al. 1999; Barvainis et al. 
2002; Cox et al. 2002) but about three times greater than those of local ULIRGs (Solomon et 
al. 1997). Assuming the most probable value for sin i=2/7i, the median ratio of gas mass to 
dynamical mass in the five galaxies is ~0.5, again 3 times greater than in ULIRGs (Downes 
& Solomon 1998), and similar to the z—2.72 Ly-break galaxy cB58 (Baker et al. 2003). Four 
of the five systems are composite AGN/starburst galaxies in a complex environment, such 
as a merger /inter acting system. The fact that the submm galaxies are complex systems is 
the more noteworthy as their redshift range is close to the peak of the merging assembly 
history of galaxy evolution. Perhaps the multiple nature of the SCUBA sources, along with 
the action of winds and outflows may explain how Chapman et al. (2003) were able to see 
strong UV line emission in sources as rich in gas and dust as the submm-population. 
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Relative to their gas reservoir, submm galaxies are very efficient emitters of radiation. 
The ratio of IR luminosity to gas mass in our five sources has a median of 380 ± 170 Lq/Mq, 
similar to high-^ QSOs (750 ± 350), HzRGs (260 ± 70; Papadopoulos et al. 2000) and local 
ULIRGs (260 ± 160; Solomon et al. 1997), but significantly larger than local LlRGs and 
more moderate luminosity starbursts (45 ± 30; Solomon et al. 1997). A young starburst 
with a Salpeter IMF between 1 and 100 Mq has ~ 10^ Lq/Mq. If most of the IR luminosity 
of our submm sources is due to star formation, their star formation efficiency must be high, 
or the IMF must be biased toward high mass stars. 

Our observations strengthen the conclusion (e.g. Genzel et al. 2003) that the brightest 
submm galaxies (Sgso ~ 2-10 mJy, corrected for lensing) have dynamical masses within the 
central few kpc that are comparable to massive, local early type galaxies (Mgas+Mgtars ~ 
m* ~ 7 X IO^^'Mq; Cole et al. 2001). Keeping in mind that our observations strictly give 
only upper limits or rough estimates of source sizes, we obtain a median dynamical mass of 
5.5x10"*^° M0 ~ 0.8 m*, for sin^^ i=7r^/4. Current semianalytic models of star formation in 
hierarchical CDM cosmogonies have difficulties accounting for the observed space density of 
such massive baryonic systems at 2; ~ 3. These models predict too few m* galaxies at that 
redshift, by about an order of magnitude, perhaps as a result of too slow baryonic cooling 
and low star formation efficiencies in the models (Genzel et al. 2003). 

We are grateful to Prof. M.Grewing for granting us the discretionary time that made 
this project possible on a short time scale. We also thank L.Tacconi for help with the data 
reduction, S.Seitz for discussions on the lensing model for SMMJ04431+0210, and I. Small, 
A.Baker and D.Lutz for thoughtful comments. We also thank O.Almaini and C.Willott 
for permission to use their optical images of N2 850.4. AWB acknowledges partial funding 
support by NSF grant AST0205937, and TRG from the EU RTN Network POE. 
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Table 1. Observed Properties for the Three SMM Sources. 









R.A. 


Decl. 


^CO 


Flux 


Line Width 


Source 


Transition 


Redshift 


(J2000.0) 


(J2000.0) 


(Jy kms~^) 


(mJy) 


(lims-i) 


SMM J04431+0210 


CO(3-2) 
CO (7-6) 
3.0 mm 
1.3 mm 


2.5094 ± .0002 


04 43 07.25 


02 10 23.3 


1.4 ± 0.2("'=) 
< 0.8 


< 0.3 
1.1 ± 0.3 


350 ± 60 


SMM J09431+4700 • • • 


CO (4-3) 
CO(9-8) 
2.8 mm 
1.3 mm 


3.3460 ± .0001 


09 43 03.74 
09 43 03.69 


47 00 15.3 
47 00 15.5 


1.1 ± O.l'"'') 
< 1 


< 0.4 
2.3 ± 0.4 


420 it 50 


SMM J 16368+4057 ■ • • 


CO (3-2) 


2.3853 ± .0014 


16 36 50.43 


40 57 34.7 


2.3 ± 0.2("'*) 




840 ±110 




CO(7-6) 


~2.383 


16 36 50.41 


40 57 34.3 


1.1 ± 0.2 








2.9 mm 










< 1 






1.3 mm 




16 36 50.40 


40 57 34.2 




2.5 ± 0.4 





Note. — The astrometric accuracy is < 0.3" (= a); Limits on Iqo and Flux are 2a, and uncertainties include statistical errors, as 
well as absolute flux errors; Line Width is full width at half maximum (FWHM); Line Velocity and 1 a error rounded to 10 kms~^; " 
No continuum subtracted; Gaussian fit; ° Frayer et al. (2003) report an upper limit of 2.5 Jy kms~^ using the OVRO interferometer. 
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Table 2. Derived Properties for the Three SMM Sources. 





Da 




D = 1'"^ 


T f a 
^ CO 






Source 


(Gpc) 




(Kpc) 


(10i°Kkins-ipc-i) 


(IO^oMq) 




SMM J04431+0210 . . . 


1.66 


4.4*^ 


1.8 


1.0 ± 0.2 


0.8 


0.3 


SMM J09431+4700 . . . 


1.53 


1.2^= 


6.2 


2.7 ± 0.3 


2.2 


1.7 


SMM J16368+4057 . . . 


1.68 




8.1 


6.9 ± 0.6 


5.5 


1.6 



Note. — Adopting a flat cosmology of VLm = 0.3, f^A = 0.7 and H = 70 kms"-*^ Mpc"'^; ° Corrected for 
the lensing magnification hl, if applicable; From Small et al. 1999; From Cowie et al. 2002; Adopting 
a conversion factor a = Mgas/i'co = 0.8 Mq (Kkms~^ pc^)~^. L'co is the apparent CO line luminosity 
corrected for the lensing magnification (Solomon et al. 1997); *= Assuming equal flux amplification and linear 
magnification; Obtained from 850/Lim flux densities (Blain et al. 2002). 
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Redshift 
2.50 2.51 2.52 

SMMJ04431 +21 
-T C0(3-2) 



3.34 



Redshift 
3.35 



3.36 



SMMJ09431+4700 
C0(4-3) 




Redshift 
2.38 2.39 

^ 1 — 



SMMJl 6368+4057 
C0(3-2) 
C0(7-6) 




1000 

LSR velocity (l<m/s) 



1000 



-500 500 

LSR velocity (km/s) 



1000 -500 500 

LSR velocity (km/s) 



1000 



Fig. 1. — CO spectra of the three SCUBA sources. The LSR velocity scale is with respect to 
the CO redshift listed in Table 1. Optical redshifts (arrows, horizontal bars are uncertainties) 
are from Frayer et al. 2003 (Hq;+[Nii], left panel), Ledlow et al. 2002 (Lya, center) and 
Small et al. 2003 (UV photospheric + Seyfert 2), respectively. The rms noise per 20 MHz 
channel is 0.7, 0.9 and 0.6 mJy for the three sources, respectively (from left to right), in 
the spectral region where the frequency settings were overlapping, and increases to about 
20 percent toward the edges of the bandpass. Overplotted on the CO 3-2 spectrum of 
SMMJ16368-I-4057 is the CO 7-6 spectrum scaled down to one-tenth of its flux density. The 
rms noise per 40 MHz channel is here 2.2 mJy. 
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C0(3-2) SMMJ04431+0210 C0(4-3) SMMJ09431 +4700 C0(3-2) SMMJ 1 6368 + 4057 




RA (J2000) RA (J2000) RA (J2000) 



Fig. 2. — Velocity integrated, natural weighted CO maps of the three SCUBA sources, 
superposed on greyscale images of the optical emission. The contours are in units of 2(T 
of the noise level, and are 0.26, 0.31 and 0.20 Jykms"^ in the three panels (left to right), 
respectively. The synthesized beams are (left to right, shown as hatched ellipses) 5.6x3.3" 
at position angle 23° (E of N), 6.6x 3.6" at 108°, and 3.3x2.6" at 79°. The three underlying 
images are in the i^-band (left panel: Frayer et al. 2003, and right panel: Small et al. 
(2003) and in the /-band (center panel: Ledlow et al. 2002). The asterisk (left image) is the 
position of the ERO N4 (uncertainty ±0.5"; see text), the filled squares (black and light) are 
the miUimeter continuum positions (center and right). The edge on spiral galaxy 2" SE of 
the CO source is the source Nl in the foreground cluster at redshift 2;=0.18. In the middle 
panel the positions of the two radio sources H6 and H7 are denoted by arrows. The stronger 
optical and radio source H6 is a narrow line Seyl galaxy at redshift 2;=3.349 (Lya). The CO 
4-3 emission from N2 850.4 remains largely unresolved. 



